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Tidal Capture of Stars by Intermediate- Mass Black Holes 



in 
o 
o 

(N 

> 
O 

(N 



> 

(N 

in 



in 
o 

Oh 
6 



H. Baumgardt 1 *, C. Hopman 2 *, S. Portegies Zwart 3,4 * and J. Makino 5 * 

1 Sternwarte, University of Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany 

2 Faculty of Physics, Weizraann Institute of Science, P.O. Box 26, Rehovot 76100, Israel 

3 Astronomical Institute "Anton Pannekoek," , University of Amsterdam, Kruislaan 403, Netherlands 

4 Section Computational Science, University of Amsterdam, Kruislaan 403, Netherlands 

5 Department of Astronomy, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 



Accepted ????. Received ?????; in original form ????? 



ABSTRACT 

Recent X-ray observations and theoretical modelling have made it plausible that some 
ultraluminous X-ray sources (ULX) are powered by intermediate-mass black holes 
(IMBHs) . N-body simulations have also shown that runaway merging of stars in dense 
star clusters is a way to form IMBHs. In the present paper we have performed TV-body 
simulations of young clusters such as MGG-11 of M82 in which IMBHs form through 
runaway merging. We took into account the effect of tidal heating of stars by the 
IMBH to study the tidal capture and disruption of stars by IMBHs. Our results show 
that the IMBHs have a high chance of capturing stars through tidal heating within a 
few core relaxation times and we find that 1/3 of all runs contain a ULX within the age 
limits of MGG-11, a result consistent with the fact that a ULX is found in this galaxy. 
Our results strengthen the case for some ULX being powered by intermediate-mass 
black holes. 
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;> , 1 INTRODUCTION 

•i-H . 

, Ultra-luminous X-ray sources (ULX) are point-like X-ray 
5_l ■ sources with isotropic X-ray luminosities in excess of L = 
10 40 erg/sec. Various theories have been proposed in the 
literature concerning the nature of ULX, like stellar-mass 
b lack ho le binaries with mild geometrical b eaming (King et 
al l 20011 . Rappaport, Podsiadlowski & Pfahl 2005|i or high- 
redshift quasars (Gutierrez & Lopez-Corredoir jj2005j) and 
it is likely that ULX are not a hom ogeneous class of ob- 
jects (Soria, Cropper & Motct i 20011) . Although luminous 
X-ray sources have been observed in different kinds of envi- 
ronments, sources with luminosities of L « 2 ■ 10 39 erg/sec 
or brighter seem to be associated with young stellar pop- 
ulations of star-for ming g alaxies (Irwin, Bregman & Athey 
I 2004 Swartz et al l 20041) . 

Since the Eddington luminosity of a star of mass M is 
given by 

M 



with luminosities exceeding L > 10 40 erg s 1 could be 
inte rmedi ate-mass black holes (IMBHs, see Miller & Col- 
be H 20041 for a review). The starburst galaxy M82 for ex- 
ample hosts a ULX with brightn ess in t he range L = (0.5 — 
1.6) - 10 41 ergs" 1 (Matsumoto et a j 200 ll. Kaaret et al l 200 ill , 
corresponding to a black hole with mass 350 — 12OOM0 if 
emitting photons at the Eddington luminosity. The case for 
an IMBH in M82 is supported by a 54mHz quasi-periodic os- 
cillation found in the X-ray flux (Strohmayer & Mushotzky 
I 2003f) and also by t he sof t X-ray spectrum of this source 
(Fiorito & Titarchul fiooll . Another argument supporting 
the connection between ULX and IMBHs are cool thermal 
emission components which have been found in the X-ray 
spectra of some ULX and which can be fitted well with ac - 
cretion disc models of IMBHs (Miller, Fabian & Mille l 20041) . 



= 1.3 x lO^ergs"' — 

Mr 



(1) 



where M is the mass of the accreting object, most low- 
luminosity ULX are probably stellar-mass black holes. How- 
ever, there is mounting evidence that the brightest ULX 
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Portegies Zwart et al. I (20041) have performed A^-body 
simulations of the dense star cluster MGG-11 whose posi- 
tion coincides with the ULX in M82. They found that run- 
away merging of massive main-sequence stars in the centre 
of this cluster leads to the formation of a massive star with 
a mass of several 1000 Mq within a few Myrs if the initial 
concentration of the cluster is larger than that of a King 
Wo = 9.0 model. The connection between an IMBH and 
the ULX was str engthened by Hopman, Portegies Zwart & 
Alexander 1(20041) who showed through analytic estimates 
that an IMBH with a mass of 1000 Mq residing in the cen- 
tre of MGG-11 has a high probability of tidally capturing a 
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passing main-sequence star or giant. The orbit of the cap- 
tured star was found to circularise even when dynamical 
perturbations by other cluster stars were taken into account 
and the system entered a Roche-lobe overflow (RLOF) phase 
once the captured star had sufficiently increased its radius 
through stellar evolution. 

P ortegies Zwart, Dewi & Maccarone If2004l) and Li ( 
2004) showed that IMBHs with mass-transferring compan- 
ions near the end of their main-sequence lifetimes are able 
to create X-ray luminosities in excess of L — 10 40 ergs _1 
over several Myrs, giving them a high chance to be observed 
as ULX. For sufficiently low-mass companions, the RLOF 
mass transfer stage might start after the host cluster has 
been dissolved by the tidal field of the parent galaxy, of- 
fering a way to e xplain isolated ULX (Hopman, Portegies 
Zwart & AlexandefjOO^). 

In this paper we present the results of detailed iV-body 
simulations of the formation and further evolution of an 
IMBH in a cluster with characteristics similar to MG G-11 i n 
the starburst Galaxy M82 (see Portegies Zwart et al l 2004) . 
Our simulations included the effects of tidal heating and the 
emission of gravitational waves of passing stars and we in- 
vestigated the chances that MGG-f I produces a ULX within 
the age limits determined by observations (7-12 Myrs). 

This paper is organised as follows: In section 2 we re- 
view the theory of tidal heating and section 3 gives analytic 
estimates of tidal inspiral. Section 4 describes the imple- 
mentation of tidal heating into our simulations. Section 5 
discusses our results for MGG-11 and section 6 summarises 
the paper and presents our conclusions. 



is given by the tidal radius (see e.g. Kochanel TT992Tl : 



2 TIDAL HEATING OF STARS NEAR A 
MASSIVE BLACK HOLE 

Stellar systems are usually well described as a system of 
point particles interacting with Newtonian gravity. In a few 
cases, however, when the stellar density is sufficiently high, 
this approximation is no longer accurate and fails to describe 
some of the important processes which occur in stellar clus- 
ters. General relativistic effects for example become impor- 
tant when stars approach each other at distances close to 
their Schwarzschild radius, which can be relevant for com- 
pact stars. Main sequence stars are much larger and close 
encounters can lead to hydrodynamical interactions, with 
stellar collisions as the extreme. This may lead to a runaway 
merger in young and dense star clus ters (P ortegies Zwart et 
al l 19991: Porte gies Zwart & McMilla Jiool Freitag, Giirkan 
fc Rasi a 20051 : see also section 5). 

Stars that pass each other at somewhat larger distances 
may still affect each others orbits by energy dissipation in 
a tidal interaction. Such an interaction may initiate tidal 
modes on both stars. The energy invested to raise the tides 
in the st ars is taken from their kinetic energies (Press & 
Teukolsk l 19771) . Two initially unbound stars may so become 
bound, or a bound eccentric binary can become bound more 
tightly. 

Stars can also have a tidal interaction with a black hole 
(BH). The length-scale for tidal forces of an IMBH of mass 
Mm on a star of mass M* and radius R* to become important 



'M.\ 1/3 
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(2) 



We now recognise several distinct regimes: (1) The 
pericenter distance is much larger than the tidal radius 
(r p 3> r t ), in which case tidal effects are negligible. (2) The 
distance of closest approach is smaller than the tidal radius 
r v < rt- In this case the entire star may be disrupted or part 
of its envelope may be stripped off to form a temporary ac- 
cretion disk around the black hole. The gas accreting from 
the star or the disc onto the black hole then causes a flare. 
The duration of the flare depends on how violent the interac- 
tion was and whether or not an accretion disc was formed: 
for tidal disruptions around super-massi ve blac k holes the 
flare might last at most a few years (Reei t 19 88). (3) In the 
intermediate regime (r p > rt), the star is tidally deformed, 
and orbital kinetic energy is transferred to the internal en- 
ergy of the star. We call this the tidal capture regime. 

For highly eccentric orbits (e > 0.9), the energy dissi- 
pation per orbit can be parametrised by 

21+2 

Ti(v), (3) 



R* 



where 



n : 



M* + Mm 



1/2 



R, 



3/2 



(4) 



Here Ti(jf) is the dimensionless tidal coupling function, 
which depends on the stellar structure and is a strongly 
decreasin g func tion of the pericenter r v (see e.g. Press & 
Teukolsk l 19771) . The orbital energy at the tidal radius usu- 
ally exceeds the binding energy of the star by several orders 
of magnitude, and some fraction of the orbital energy is dis- 
sipated at every new pericenter passage. As a result the star 
becomes very hot and expands. Both effects make it more 
luminous. The stellar luminosity can become much larger 
than that achieved by nuclear burning, and is of order 

Lt = — • (5) 

Here P is the orbital period of the binary system. 

It is unclear where a star stores the excess tidal energy. 
In the two most extreme cases, the star can store the ex- 
cess energy in the surface layers, or in the bulk of the stellar 
material. The first case leads to high surface temperatures 
without appreciable expansion of the star (McMillan, Der- 
mott & Taar rl 19871) . whereas the latter causes the star to 
remain cold but exp and dr amatically (Podsiadlowsk F"l996l) . 
Alexander & Morri j 20031 dubbed the two possibilities as 
"hot" and "cold" squeezars and argued that these tidally 
excited stars may be observable in the center of the Milky- 
Way Galaxy. 

The long-term response of the star to tidal heating de- 
termines whether the star survives the encounter. To be able 
to survive the tidal circularisation near r t , a star on an ini- 
tially wide orbit has to dissipate 



Et 



M. \ 2/3 GMl 



R* 



(6) 



which generally is much larger than the binding energy of 
the star. As a consequence, a star can only survive a strong 
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tidal interaction if it cools efficiently. The most obvious way 
this can be achieved is by radiation, i.e: a "hot squeezar". 

The structure of a "hot squeezar" is unaffected by the 
encounter, and as a result the cooling of the star is limited by 
the Eddington luminosity. The limiting luminosity can now 
be used to compute a lower limit on the pericenter for which 
the star is able to survive an encounter. The further orbital 
evolution can be described analytically and the time-scale 
for com plete c ircularisation of the orbit is then (Alexander 
& Morri i 2003ft : 



t (r p ,a) = 



2irMWGM.a 

AE t (r p ) 



(7) 



If this time-scale is sufficiently small, the orbit of the 
star is not changed significantly by scattering with field stars 
during inspiral (section l3.3H . In this case the star circularises 
near the tidal radius of the IMBH. The subsequent evolution 
of the IMBH with main-sequence star binary is discussed in 
section 15.51 . 



3 ANALYTICAL ESTIMATES OF TIDAL 

INSPIRAL IN PRESENCE OF SCATTERING 

If a star has a orbit in the tidal capture regime (r p > rt) af- 
ter a first encounter with an IMBH, it experiences repeated 
tidal interactions through which the orbit shrinks and and 
becomes circularised. The time-scale for circularisation is 
much larger than the orbital period P of the star (to ^> P). 
If the orbit is perturbed by another star during the inspi- 
ral process, the inspiraling star may either be scattered to a 
wider orbit or into a tighter orbit. In the first case tidal heat- 
ing will become less effective and the inspiral process slows 
down or stops completely, making the star more vulnera- 
ble for subsequent dynamical interactions. In the latter case 
the inspiral process is either accelerated or, if the pericenter 
distance becomes smaller than the tidal radius of the star, 
the star is tidally disrupted by the IMBH and the inspiral 
process ends. 



3.1 Assumptions 

We present an approach which is quite similar for tidal dis- 
ruptions and tidal inspirals. To obtain analytical estimates 
for the two processes, we make several assumptions that are 
generally made for th e treat ment of tidal disruptions (e.g. 
Lightm an and Shapir d 19771: Fran k & Ree J 1976t Cohn & 
Kulsr u Jl97l Syer & Ulm d 1999| M agorr ian fc T remaine 
I 19991: Miralda- Escud e & Goul d 20001: Fr eitat 2001al Alexan- 
der & Hopmar pOOl Wan g fc M erri d 2004ft "and i nspiral 
processes fe.g. Hils fc Bendel- 19951 : Sigurdsson & Rees 1997; 
Ivano f200a Freita l 2001 bl 20031: Alexan der fc H opma £ 2003t 
Hopman, Portegies Zwart & Alexande J 20041 : Hopman & 
Alexande F2005Tl . 

Our simple analytical model captures some of the phys- 
ical mechanisms of the cluster and gives approximately the 
correct estimates for the probability that a star is captured 
by the IMBH. Since many of the assumptions listed here 
have been made by numerous previous papers, it is of im- 
portance to trace to what extend the analytical model agrees 
with the simulations, which has less number of simplifying 
assumptions. 



Here we list our assumptions and discuss them. 

(i) The stellar distribution function is well approximated 
by a stellar cusp, i.e. 



*(r) = 



(3 - a)N a 



(8) 



where a is approximately constant withi n the radius of in- 
fluence r a of the IMBH (Bahcall & Wol f 197fill 1977ft . and 
N a is the number of stars within r a . 

This assumption is shown t o be satisfied by earlier N-body 
simulation s (Bau mgardt et al l 2004 J 2004b: Preto, Mcrritt 
& Spurzem 2004), and also in our current simulations. The 
slope of the cusp in our simulations is a ~ 1.5. A stellar cusp 
with a « 1.5 has also been observed by s tar co unts near 
the MBH in our Galactic Centre fAlexande f 1999ft . We will 
assume a = 3/2 in the following analysis. This simplifies the 
expressions somewhat, because in that case the relaxation 
time is constant (see eq. 1121 1. 

Baumgardt, Makino & Ebisuzaki: (20 04ah have shown that 
for small-mass black holes, containing less than a few percent 
of the total cluster mass, the radius of influence of the black 
hole is given by 



GM. 



(9) 



where a is the velocity dispersion in the cluster core. 

(ii) Within r a the orbits of the stars are Keplerian to good 
approximation. This is a good assumption for our purposes, 
and it makes many of the expressions more transparent. 

(iii) The relaxation time within the cusp is much larger 
than the orbital time. Stars exchange energy and angular 
momentum only through small angle two-body interactions. 

(iv) Within the radius of influence the stellar population 
can be approximated by a single mass population. 

This assumption is made for simplicity, and much of the 
discrepancy between the analytical model and the simula- 
tions stems from this assumption. In a young stellar cluster 
there is a wide range of masses, with a few stars which have 
masses much larger than the mean stellar mass (M*). As 
a result the most massive stars sink to the IMBH, and the 
masses within the radius of influence vary strongly. This has 
considerable consequences for the dynamical behaviour close 
to the IMBH which the simple analytical model fails to de- 
scribe. We will comment on this when we discuss the results 
of our simulations. 

(v) Stars reach the IMBH by diffusion of angular momen- 
tum rather than energy. 

This assumpt ion is a ccurate for a single mass distribution 
(Bahcall & Wol f 1976ft . When a few very massive stars are 
present in the cluster, however, dynamical friction operates 
on a time-scale ia.f. ~ ((M+) / M*)t r , much shorter than the 
relaxation time t r . The dynamical mechanism which drives 
stars to inspiral orbits in a young stellar cluster therefore 
typically operates in two phases. First there is an "energy 
phase", during which a massive star loses energy to field 
stars and sinks to the centre; this is a purely elastic phase 
during which there is no tidal heating. Then follows an "an- 
gular momentum phase", during which more interactions 
with cluster stars change the angular momentum of the mas- 
sive star, until it finally reaches an orbit with pericenter close 
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enough to the tidal radius, and tidal heating becomes effi- 
cient. 

As a result of our simplifying assumptions, the following 
analysis does not describe the radial distribution function 
(DF) of the stars as a function of their mass. But once this 
DF is given, it does provide a treatment of the interplay be- 
tween scattering and dissipation which is in good agreement 
with the simulations. 



3.2 Tidal disruption 

Tidal disruption of stars by an IMBH occurs when a star 
on an eccentric orbit has angular momentum smaller than 
the angular momentum at the "loss cone" (J; c ). For a main- 
sequence star orbiting an IMBH, the loss-cone is defined as: 



Jic = V2GM.r t . 



(10) 



The disruption rate and distribution of stars in an- 
gular momentum were studied by Lightman & Shapiro ( 
1977). Here the average changes in momentum are much 
smaller than Ji c (diffusive regime). For this case the time- 
scale for changes of the p ericent er by (small angle) scattering 
is (Alexander & Hopmar l 20031) 

t p (r p , a) — —t r . (11) 
a 

Here t r is the relaxation time, which we can estimate with 

P 



tr=A A [ — 



N(< r) 



(12) 



Here Aa is a constant which includes the Coulomb loga- 
rithm, and 7V(< r) tx r s ~ a is the number of stars enclosed 
within r. For a = 3/2 the relaxation time is independent of 
distance from the IMBH. 

The diffusive regime is then defined to be the region 
where a < r cr i t ; r- cr it is the semi-major axis for which 
-P(rcrit) = t p (r t , r cr it), 



2?r 



VGM.r t t, 



2/5 



(13) 



The angular momentum of stars with a < r cr i t changes by 
an amount smaller than J; c per orbit, while the change in 
angular momentum per orbit is larger than Ji c if a > r cr it. 
In the latter case the velocity distribution can be isotropic, 
since it is not affected by the presence of the loss-cone. This 
is not true for the diffusive regime. 

The distribution function and merger rate in the diffu- 
sive regime are computed by solving the appropriate Fokker- 
Planck equation. The distribution function should become 
isotropic for J 3> J; c and vanish at the loss-cone. The steady 
state solution gives rise to a flow of stars through the loss- 
cone, which is a consequence of the presence of a mass sink 
at Ji c . This is reflected in the boundary condition that the 
distribution function vanishes at the loss-cone. 

The disruption rate in the diffusive regime is given by 



daN iso (a) (r clit /r a ) 3/2 N a 



i t r ln(J m / Jic) t r \n(J m / Jic) 
(e.g. Lightman & Shaping 19771: Hopman & Alexande F^OO. 



(14) 



Here N(a)da oc a 2 ~ a da is the number of stars with semi- 
major axes in the range (o, a + da). In eq. (1141 J m is to be 
evaluated at 7- cr it. 

Equation II H is interpreted as a slow diffusion of the 
angular-momentum of the stars. The time-scale for the an- 
gular momentum to change by a factor ~ J m = \J GM,a is 
of the order of the relaxation time t r . In a time t r a signifi- 
cant fraction of the stars within 7- cr it are therefore disrupted 
by the IMBH. The disruption rate is suppressed by a log- 
arithmic factor which reflects the presence of the loss-cone. 
This leads to a dilution of stars with small angular momenta 
in the diffusive regime, as can be seen from the s olution of 
the Fokker-Planck equation (Hopman & Alexander 2005). 

For the tidal disruption rate, there is also a contribution 
from stars for which the change in angular momentum per 
orbit is larger than J; c (i.e., the regime a > r cr i t ). In that 
case the disruption rate T(< a) for stars with semi-major 
axis smaller than a is estimated to be given by the fraction 
of stars in the loss-cone, divided by the period, 



I* ick (< a) 



N(< a)r t 
aP(a) 



(15) 



For tidal capture (see below) there is no such con- 
tribution, and we do not further consider stars in this 
"kick" reg ime h ere; for a discussio n see, e.g., Lightman 
fc Sh apiro [7l977h : Cohn & K ulsrud ffi97gh : Syer & Ulmer 
Ifl999h : Magorrian & Tremaine f?1999T) . 



3.3 Tidal capture 

Due to relaxation in angular momentum, any star will even- 
tually be accreted by the IMBH, either because of direct 
tidal disruption, or due to tidal heating and inspiral. If 
the number of stars in (a, a + da) which spiral in is given 
by Ni(a)da, and the number of stars which are directly 
disrupted by Nd(a)da, then the probability that the star 
reaches the IMBH by inspiral is given by 



S(a) = 



Ni(a) 



N t (a) + N d (a) 



(16) 



Hopman & Alexander If2005h performed Monte Carlo 
simulations to show that stars with high energies (small 
semi-major axes) are unlikely to enter the IMBH without 
first spiralling in, while stars with low energies (large semi- 
major axes) are tidally disrupted without experiencing a 
tidal inspiral phase. 

The rate of inspiral is then given by 



r, = 



1 daN(a)S(a) 

t r \n(Jm/Jlc) 



(17) 



Here the function S(a) to accounts for the fact that capture 
without premature disruption is improbable for a > a c . 

In order for the star to survive, tidal heating should hap- 
pen at a sub-Eddington rate (see previous section), yielding 
a minimal value for the pericenter r p , m j n ~ 2r t . Inspiral 
without premature scattering is only feasible provided that 
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to(r p ,a) < t p (r p ,a). As a consequence stars can only spiral 
in if initially their semi-major axis is smaller than 



SAEt (^*p,min)^*p, miliar 



1 2/3 



(18) 

Note that to > P and thus a max < r cl it, which implies that 
for tidal capture the relevant regime is the diffusive regime, 
and that the tidal capture rate is smaller than the rate for 
tidal disruption, in spite of the fact that tidal heating is 
efficient for r p < 3rt, while tidal disruption happens only 
for r p <r t . 

The function S(a) may be approximated by a Heavi- 
side step- functio n 8(a — a c ), where a c > a max (Hopman & 
Alexande ilooil . This approximate behaviour is confirmed 
by our simulations (see JSJ. The rate at which stars undergo 
inspiral due to tidal heating and survive is then given by 



r, 



: daTV(q) 

t r ln( J m j Jl c 



(amax/r a ) 3/2 jV a 

t r ln(J m /Ji c ) 



(19) 



(20) 



Here the radial stellar distribution function is given by Eq.|H] 
with a — 3/2. 

Comparing Eas. l20l and ll4l implies that the rate at which 
stars in the diffusive regime are tidally disrupted is larger 
than the tidal capture rate by 



r, 



^"crit 



3/2 



(21) 



The actual tidal disruption rate is even higher, because 
stars with semi-major axes a > r cr i t can also be disrupted. 



4 DESCRIPTION OF THE TV-BODY RUNS 

We simulated the evolution of MGG-11 through TV-body 
simulations of star clusters containing TV = 131,072 (128K) 
stars u sing A arseth's collisional TV-body code NBODY4 
fAarset fl999T) on the G RAPE6 computers of Tokyo Univer- 
sity (M akino et all_2003). Following Portegies Zwart et al. ( 
2004), our clusters are initially King models with Wo = 9 
and half-mass radius — 1.3 pc. The initial mass spectrum 
of the cluster stars was a Salpeter mass-function between 
I.OMq ^ m ^ 1OOM0. These cluster models have a pro- 
jected half-mass radius, mass-to-light ratio and total cluster 
mass after 12 Myrs that is consistent w ith the properties of 
MGG-11 as observed by McCrady et al. lf2003l) . The clusters 
are also concentrated enough to form IMBHs via runaway 
merging of st ars. St ellar evolution was modelled according to 
Hurley et al. [(200 J) and individual runs were calculated for 
T = 12 Myrs. Table 1 gives an overview of the simulations 
performed. 

Stars were merged if their separation became smaller 
than the sum of their radii and the mass of the merger prod- 
uct was set equal to the sum of the masses of the stars, i.e. 
we assumed no disruption of stars and no mass loss during 
a collision. Since typical velocities of stars in our clusters 
are much smaller than the escape velocities from stellar sur- 
faces, this should be a valid assumption. The mass of tidally 



disrupted stars was added to the mass of the central black 
hole. 

In total we performed 12 runs. In six of the runs we 
used a relation for the radius of the runaway star d erived 
from a fit to the results of Ishii, Ueno & Kato lfl999l) given 
by 



1 R * 1 K i M * 
log i^ = L51og Me 



1.85. 



(22) 



for stars with 120M Q < Af* < 1000M©. Radii of more 
massive stars were set equal to the size of 1O 3 M0 stars. 
These runs created fairly massive IMBHs with masses be- 
tween 3OOOM0 < m < 4000Mq. In another set of runs, we 
used t he mass-radius relation found by Bond, Arnett & Carr 



Ifl984h : 



log-gl = 0.47 log -^ + 0.20. 

ft© M e 



(23) 



While derived for stars with masses larger than 10 4 Mq, 
this relation gives radii in good agreement with the ZAMS 
radii from Hurley et al. for stars with masses around 10 2 Mq. 
These runs lead to runaway stars with masses around 
1000-2000M Q , roughly equal to the estimated mass of M82 
ULX-1. We did not apply stellar mass- loss for the runaway 
stars and transformed them into IMBHs at T = 3 Myrs. 



4.1 Implementation of tidal heating into the 
TV-body simulations 

We have implemented orbital energy loss by tidal interac- 
tions between a star and the IMBH in the TV-body simula- 
tions using the prescription of Portegies Zwart & Meinen ( 
I 19931 . henceforth PZM). They present fitting formula to the 
tidal energy loss as a function of the masses, radii and poly- 
tropic indices of passing stars. The dimensionless functions 
Ti(rf) in eq. @ for the tidal energy were obtai ned fr om a 
polynomia l fit to calculations by Lee & Ostrikcr (1986) and 
Ray et al. Ifl987j) for stars of polytropic indices n = 1.5,2 
and 3. PZM found that the / = 2 and 1 = 3 terms give a 
99% contribution to the dissipated energy and therefore did 
not take higher harmonics into account. 

In our implementation, we assumed that main sequence 
stars with masses m < O.4M0 and giants have polytropic 
index n = 1.5, while stars with mass O.4M0 < m < 1.25M0 
have polytropic index n — 2. More massive main-sequence 
stars were assumed to have n — 3. It was not necessary to 
apply tidal heating to white dwarfs or neutron stars since 
these stars were not present in the runs because the studied 
clusters are too young. 

Eq. 3 is only valid if the orbital eccentricity is close to 1. 
For orbits of smaller eccentricity the amount of tidal heating 
is smaller since the star feels the influence of the IMBH also 
at apocenter. In order to reflect this we multiply the heating 
energy AEt by the eccentricity e of the orbit. In our current 
implementation, the total angular momentum is conserved 
and an unperturbed IMBH-star system will always circu- 
larise at twice the initial pericenter distance r P i. 

Due to the structure of NBODY4, tidal heating had to 
be included in three different places in the code. In the main 
integrator, we checked the distance of each particle to its 
nearest neighbour each time it was advanced in time. This 
is facilitated by the GRAPE hardware which provides the 
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Table 1. Details of the performed N-body runs. Shown are the mass-radius relation assumed for the runaway stars (Ishii et al. or Bond 
et al.), the mass M. j of the IMBH at time of its formation (T=3 Myrs), the mass M, f of the IMBH at the end of the run (T=12 Myrs) 
and the number Njjis of stars tidally disrupted by the IMBH. Final columns give orbital parameters for tidal heating events: Starting 
time Ti nsp of the inspiral, duration of inspiral ATj nsp until the orbit is circularised, mass M* and radius ij* of the star at time of inspiral, 
tidal radius Rt of the star, initial pericentre distance of the orbit r p j, ratio of pericentre distance to the tidal radius of the star r p ;/r t , 
semi-major axis a and eccentricity e of the initial orbit and semi-major axis r Pi f after the orbit circularised. 
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index of the nearest neighbour for each particle at almost 
no extra computational cost. The closest encounter distance 
was calculated once the relative velocity v re i — ^ij^ij 
between two particles switched sign from negative to pos- 
itive and we applied tidal heating to all stars that passed 
within 5r t of a massive black hole by decreasing the kinetic 



energy of the combined star-black hole system by an amount 
AE t . Energy was extracted by changing the velocity and po- 
sition of the black hole and the star when the system was 
at pericenter. The energy was changed according to eq. 3 
such that the total linear and angular momentum remains 
constant. 
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Figure 1. Density profile of run 6 at T = 3.15 Myrs, just before 
an inspiral event happens. Between 10 pc < r < 0.05 pc, the 
density profile of the cluster follows a power-law cusp with a = 
1.5. The innermost data point is due to one star which experiences 
an inspiral at T = 3.15 Myrs. 



Since in NBODY4 the motion of close particle pairs 
is foll owed by KS regularisation (Kustaanheimo & Stiefel 
I 19651) . we also had to add tidal heating to the motion of 
KS binaries. For unperturbed pairs this was done each time 
their motion was integrated by changing the semi-major axis 
and eccentricity of the binary according to the accumulated 
amount of tidal energy. For perturbed binaries, tidal heat- 
ing was applied after each pericenter passage since at this 
point the perturbation of the motion of close neighbours is 
smallest and it is easiest to correct the total energy of the 
cluster. 

Energy dissipation by grav itational radiation was imple- 
mented following Peters (1964) formulae in a similar fashion 
as the implementation of tidal heating. 



5 RESULTS 

5.1 Central density and mass segregation 

We begin by discussing the density profiles of the star clus- 
ters after an IMBH has formed in the centre. Fig. 1 depicts 
the density profile of run 6 from Table 1 at T = 3.15 Myrs 
just before a tidal capture event occurs. The cluster centre 
is assumed to coincide with the position of the IMBH. It 
can be seen that the cluster follows a power-law mass den- 
sity profile p ~ r~ a with slope a — 1.50 inside r = 0.05 
pc. This is close to the va lue of a = 1.55 found by Baum- 
gardt, Makino & Ebisuzaki[?2bo4b) for IMBHs in multi-mass 
clusters. Inspection of the calculations shows that a cusp is 
already present by the time the runaway star has collapsed 
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Figure 2. Average and maximum masses of stars in the cluster 
as a function of time. The masses of all stars decrease as a result 
of stellar evolution. Due to mass segregation, stars in the cusp 
around the IMBH are more massive than the average cluster stars 
and the star closest to the IMBH is always one of the most massive 
stars in the cluster. 



to a black hole (at about 3 Myrs), and as a result already 
influences the late stages of runaway merging. One reason 
for the quick formation of cusps is the large initial mass 
spectrum of our models which reduces the relaxation time. 
Fig. 1 also justifies our assumptions about the presence of 
a cusp around the IMBH in §3. Due to the finite number of 
cluster stars, the central cusp runs out of stars at r ~ 10~ 3 
pc, corresponding to r = 10 5 Rq. In the present example, 
one star was pushed to r = 10~ 4 pc due to dynamical inter- 
actions between the inner cusp stars. This star experiences 
an inspiral event at T = 3.15 Myrs. 

Fig. |21 shows average and maximum masses of stars in 
the cluster and near the IMBH over the course of a simu- 
lation. Values shown are averaged over several simulations. 
Masses of high-mass stars decrease as a result of stellar evo- 
lution since stars lose mass in winds and the most massive 
stars are constantly removed and turned into compact rem- 
nants. As a result, the mass of the most massive cluster star, 
which is near 80 Mq at T = 3 Myrs, drops to 15 Mq at the 
end of the simulation. Due to the large number of low-mass 
stars, the average mass of cluster stars stays nearly constant 
within the first 10 Myrs at (m) ~ 3M®. Due to mass seg- 
regation, the mean mass of stars in the cusp is higher than 
the mean mass of all cluster stars. Also within the cusp, 
massive stars sink close to the IMBH so the stars near the 
IMBH have masses significantly higher than the mean mass. 
This shows that it is important to take mass segregation into 
account when determing the inspiral rates. 

5.2 Tidal capture 

Several trends can be deduced from Table 1: 



8 H. Baumgardt, C. Hopman and S. Portegies Zwart, J. Makino 



• Successful inspiral events occur in all runs with on av- 
erage ~ 4 inspiral events per run during the first 9 Myrs 
after IMBH formation. Hence there is a high chance that an 
IMBH in a young star cluster like MGG-11 will undergo a 
RLOF phase at least once during the lifetime of the cluster. 

• Since the masses of stars in the inner cusp are relatively 
large (see Fig. |2J , most inspiral events also involve massive 
stars. The average mass of captured stars is (m) w 25Mq, 
which exceeds the average mass of stars in the cusp. This 
is caused by the larger cross-section for a tidal interaction 
with the IMBH of the massive stars. 

• In agreement with the theoretical estimates of §3, suc- 
cessful tidal inspiral occurs only for stars in tightly bound 
orbits around the IMBH with semi-major axis a < 10 5 -Rq, 
the exception being giant stars with much larger stellar radii. 
This corresponds roughly to the inner end of the cusp. Stars 
on orbits with larger semi-major also experience tidal inter- 
actions with the IMBH (see Fig. |7J , but are scattered away 
by cusp stars before their orbits can circularise. 

• Although the cross-section for tidal heating is larger 
than the cross-section for tidal disruption, it is much more 
likely to tidally disrupt stars, than to bring a star into a close 
orbit around the IMBH via tidal heating. Averaged over all 
runs, we find that the ratio of tidal disruptions to inspirals 
is N d /Ni w 7. 

• Successful inspiral events occur predominantly early on 
in the cluster evolution: Between 5 < T < 12 Myrs only 
24 inspiral events occur compared to 26 events between 
3 < T < 5 Myrs. The decay in the number of tidal inspi- 
ral events is mainly caused by the expansion of the cluster, 
which is driven by mass loss from stellar evolution. In addi- 
tion, the most massive cluster stars leave the main sequence 
and increase their radii at T w 3.3 Myrs, which explains the 
large number of inspirals around this time involving stars 
with masses m > 25M©. 

• The final pericenter distances after the orbit circu- 
larised, are within a factor of three of the tidal radii since 
the timescale for tidal circularisation is small only for small 
r p /r t . 

As an example for an inspiral event, Fig. [3] depicts the 
evolution of semi-major axis a and pericenter distance r p of 
a captured m — 17.6 Mq star with time from one of our 
runs (run 3). The time when tidal heating shrinks the orbit 
is shown by a solid line. Between 5.3 and 6.1 Myrs the star 
segregates towards the IMBH as a result of its higher than 
average mass and the semi-major axis of the orbit shrinks 
to a few 10 4 Rq. Tidal heating occurs at T = 6.13 Myrs, 
when the pericenter distance is scattered to a distance only 
slightly larger than the tidal radius. During the inspiral the 
semi-major axis decreases from more than 10 4 R© to about 
120 Rq, and the orbit circularises at r p f/Rt = 2.5. After 
the tidal circularisation, the semi-major axis remains un- 
changed except for perturbations due to passing stars, which 
increase the eccentricity. Upon each time a small eccentric- 
ity is induced, tidal effects will directly start to circularise 
the orbit again, and so result in a slight decay of the semi- 
major axis. A strong orbital perturbation causes the star 
to fill its Roche-lobe at T ~ 10.2 Myr and it is destroyed 
by the IMBH. Our present TV-body code does not allow to 
treat RLOF phases, a star whose stellar radius expands or 
whose orbital radius shrinks to the point where the stellar 
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Figure 3. Semi-major axis and pericenter distance of a tidally 
heated star before and after inspiral due to tidal heating. The 
main tidal heating event happens at T = 6.13 Myrs at which 
time the star circularises. Later, the semi-major axis changes due 
to external perturbations while the tidal radius increases as a 
result of stellar evolution and IMBH growth. 

radius exceeds the tidal radius will simply be disrupted. We 
will present a discussion about the further evolution of such 
systems based on simple binary-stellar evolution estimates 
in chapter 5.5. 

5.3 Super-Eddington heating rates 

The inspiral process may lead to complete tidal circulari- 
sation, but also to the destruction of the star. The latter 
happens if 1) the star passes the black hole's event horizon 
or 2) if the rate of energy dissipation in the tidal encounter 
exceeds the Eddington luminosity. We will now estimate the 
latter fraction. For each inspiral event we calculated the 
rate at which energy was pumped into the star due to tidal 
heating when the eccentricity of the orbit was e = 0.9. At 
this time the amount of energy dissipation per time AE/P 
reaches a maximum since a and P decrease with e while for 
smaller eccentricities the stars feel the influence of the black 
hole also at apocenter and AE decreases. 

Fig. [I] compares the energy dissipation rates at e = 
0.9 with the Eddington luminosities of the stars. Super- 
Eddington heating rates occur in roughly half the cases, 
but since Fig. [I] shows the maximum heating rates during 
the inspiral which occur only for a few orbits, cases with 
AEi„s P > LEdd do not necessarily end up in the destruc- 
tion of the star. Stars might for example lose their outer 
envelopes and shrink, thereby becoming less susceptible to 
tidal perturbations, so that the central parts of the stars 
survive the inspiral. In addition, the extra energy pumped 
into the star must be larger than the potential energy of the 
star in order to disrupt it, so the super-Eddington heating 
rates have to be maintained over many orbits. In order to 
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Figure 4. Energy dissipation rates AEi nap when the orbital ec- 
centricities have reached e = 0.9 compared with the Eddington 
luminosities L^dd of the stars. Nearly half of all stars experi- 
ence sub- Eddington heating rates during the inspiral process and 
survive the inspiral. Moderately super-Eddington stars might also 
survive the inspiral since the super-Eddington heating rates occur 
for only short intervals and the extra energy pumped into the star 
is not sufficient to disperse it completely. If the amount of tidal 
heating is too large, the stars expand and are tidally disrupted 
by the IMBH. 



answer the question which stars survive the inspiral, more 
detailed simulations treating the response of the star to the 
tidal heating would be required. These are beyond the scope 
of the present paper. 

5.4 Comparison with analytic estimates 

We compare the numerical results with the analytical es- 
timates from for the following rather typical example: 
M. = 3 x 1O 3 M , M* = 30M Q , = WR Q , t rc = 3 Myrs, 
r a = 0.05 pc, N a = 300, and r p /r t = 2.3. With eqs. 
(14), (19) and (21) we then find that over a time inter- 
val of Tf =9 Myr the number of tidal disruptions is 
Nd = TdT/ = 8.6, while the number of tidal captures is 
smaller by a factor (r cri ,/a milI ) 3/2 ~ 0.08, Ni = YiT f ~ 0.69, 
i.e: Ni/Nd — 0.08. The maximum distance from which stars 
can originate is a max ~ 3.4 x lO 4 i?0. By using eq. 11311 . we 
also find that the critical radius inside which the loss cone 
should become empty is r cr it = 2 ■ 10 5 i?Q, slightly larger 
than the inner end of the cusp. 

From our simulations we find that Nd ~ 27, Ni = 4.1 
and Ni/Nd = 0.15. Using only the inspirals which have sub- 
Eddington heating rates gives Ni = 1.58 and Ni/Nd = 0.06 
Hence the results for the analytical inspiral and disruption 
rates Ni and Nd are within a factor of 3 to those found in 
the simulations while the ratio Ni/Nd is correctly predicted. 
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Figure 5. Initial semi-major axis distribution of tidally dis- 
rupted stars (solid lines) and combined distribution of tidally dis- 
rupted and captured stars (dotted). Shown arc differential rates 
N(a) = (N(> a)—N(> a + Aa))/Aa. The combined distribution 
has a maximum near the critical radius r cr it where the loss-cone 
becomes empty, in agreement with theoretical predictions (Light- 
man & Shapiro 1977). The combined distribution can be fitted 
by two power-laws. 



Although tidal disruptions are not the main focus of 
the present paper, it is interesting to compare them with 
analytic estimates since the physical processes responsible 
for inspirals and disruptions are the same. The distribution 
of semi-major axis of tidally d isrupte d stars has been calcu- 
lated by Lightman & Shapiro [719771) . They found that the 
peak in semi-major axis distribution should occur near r cr it 
(sec cq. 1131 ). Inside r < r cr i t , the loss cone is empty, so stars 
have to be scattered into the loss cone by two-body relax- 
ation. In this case the number of disrupted stars should be 
proportional to Nd(< r) ~ N(< r)/t r which leads to a scal- 
ing Nd(< a) ~ a 3//2 for the number of encounters with semi- 
major axis smaller than a (see eq. 1141 ). For the correspond- 
ing differential function we find therefore Nd(a) ~ a 1 ' 2 . At 
larger distance r > r cr it the loss cone is full ( "kick" regime) ; 
the number of disrupted stars is then given by eq. |(T^J. This 
leads to a scaling Nd(> a) ~ a~ , or, for the differential 
function Nd(a) ~ a~ 2 . 

Fig. shows the number of circularised and tidally dis- 
rupted stars as a function of the initial semi-major axis. The 
maximum occurs at a pca k — 2 • 10 4 Rq, somewhat smaller 
than our analytical estimate for the critical radius rcrit - 

For radii smaller than a pGa k, the combined distribution 
of tidally heated and disrupted stars can be fitted by a rela- 
tion N(a) ~ a 1,5 . This is somewhat steeper than predicted 
by Lightman & Shapiro lfl977h . however the uncertainties in 
the observed slope are quite large due to the small number 
of events at small radii. If real, the reason for the discrep- 
ancy with Lightman & Shapiro could be the capture of stars 
by the IMBH followed by three body encounters with other 
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Figure 6. Fraction of tidally captured stars compared to all stars 
ending in or near the IMBH as a function of the initial semi-major 
axis a of the stars before tidal capture/disruption (see S(a) in eq. 
[ED)- For small semi-major axis tidal capture dominates while 
for larger radii tidal disruption becomes the dominant process 
for feeding the central black hole since stars cannot spiral in any 
more before being scattered away by other stars. 



Figure 7. Semi-major axis of tidally circularised stars (dots) 
compared to stars for which tidal heating does not lead to cir- 
cularisation (crosses). Successful tidal inspiral is only possible for 
stars on orbits with small semi-major axis. The threshold between 
the two cases is in agreement with the prediction of this paper 
and of Hopman et al. (2004) for the maximum inspiral distance 



stars in the cusp. An inspection of the output data shows 
that these processes contribute considerably to the tidal de- 
struction rate but are not taken into account by our analytic 
estimates. For radii larger than a pea k, we obtain a slow er de- 
crease than predicted by Lightman & Shapiro 1(19771) . The 
reason could be that the stars in the cusp drive the Brown- 
ian motion of the IMBH. This influences the determination 
of the semi-major axis of in-falling stars in case of large 
semi-major axis. 

Hopman & Alexander I (2005|) found that if stars move 
slowly through phase space, the number of tidally captured 
stars should be much larger than the number of tidally dis- 
rupted stars for small radii, Ni(a) S> Nd(a); while for large 
radii two-body relaxation becomes important and prevents 
tidal inspiral so Nd(a) Ni(a). In Fig. |S|we compare the 
number of stars that are tidally heated with all stars that 
strongly interact with the IMBH. It can be seen that we 
indeed find the predicted behaviour for S(a). 

In §3.3 we demonstrated that the semi-major axis of a 
stellar orbit has to be smaller than a max in order for the or- 
bit of the star to be successfully circularised. Otherwise the 
star is scattered away from its orbit due to encounters with 
other stars. Fig.Qshows the initial semi-major axis of stars 
which are tidally heated by the IMBH. Stars for which the 
inspiral continues all the way until the orbit is circularised 
are shown by filled circles while stars for which the heating 
process is terminated since they are either scattered away or 
experience super-Eddington heating rates during the inspi- 
ral process are shown with crosses and triangles respectively. 

Tidal capture is only successful if the initial semi-major 
axis of the star a < 2 ■ lO 5 i? . The average inspiral distance 



agrees quite well with the predicted value of a max . Stars cir- 
cularising due to tidal heating fall into two categories: most 
of them initially moved on orbits with e « 0.99 and come 
from the inner cusp, while few came from relatively low- 
eccentricity orbits with semi-major axis of a few 1000-Rq. 
Most of the stars originating from low-a orbits experienced 
dynamical interactions with previously inspiraled stars. 

5.5 Further Evolution of captured Stars 

In order to study the mass accretion rate and X-ray luminos- 
ity of the IMBHs, we assume that only the innermost stars 
are close enough to transfer a significant amount of mass to 
the IMBHs and u se the model discussed by Kording, Falckc 
& Markoff[(2002). They assume that X-rays are created in 
accretion discs and jets and argue that a black hole with a 
nearby companion is in the hard state if M > M C rit in which 
case the disc luminosity is given by Lx = eMc 2 . At lower 
accretion rates Lx = eM 2 c 2 /M CT i t . Following their paper, 
we adopt e — 0.1 and Af cr it = O.1M0/Myr. The mass-loss 
rates of the donor stars are direc tly tak en from the stellar- 
evolution routines of Hurley et al. (2000). Stellar winds have 
typical velocities of about 200 km/sec for LBVs and 1000 
km/sec for Wolf-Rayet stars (Paumard et al. 2001), so de- 
pending on the orbital separation of the innermost star and 
the IMBH, gas lost by the star can also escape from the sys- 
tem instead of ending up on the IMBH. A detailed study of 
gas accretion would require SPH c alculations like the ones 
presented by Cuadra et al. [(2005) , which are beyond the 
scope of this paper. In order to account for the loss of gas, 
we have plotted times when the escape velocity of the star 
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Figure 8. Mass accretion rate and X-ray luminosity of the IMBH 
of run 10. Crosses denote times when the distance of the star 
closest to the IMBH is so large that most of the mass lost from 
the star escapes from the star-IMBH system. Between 3 and 12 
Myrs, several close companions to the IMBH formed through tidal 
heating which transfer mass to the IMBH. For two events, the X- 
ray luminosity of the IMBH reaches Lx = 10 40 erg/sec while 
the stars are on the main-sequence, enough to explain the X-ray 
luminosity of M82 ULX-1 and also within the right age range of 
MGG-11. 

is smaller than 200 km/sec with crosses. In these cases the 
values for Lx should be regarded as upper limits. 

Figs.|H|and|n|depict the X-ray luminosities of the IMBHs 
as a function of time for runs 7 and 10. For both clusters, the 
distance and nature of the innermost star changes rapidly 
during the simulation as stars near the IMBH are turned into 
compact remnants by stellar evolution and then scattered 
out of the cusp by other stars or are tidally disrupted by the 
IMBH. Run 10 experiences 4 prolonged ULX phases between 
3 and 12 Myrs, each time due to massive main-sequence stars 
brought close to the IMBH as a result of tidal heating. The 
last two phases are within the right age limits for MGG-11 
as determined by McCrady et al. and produce X-ray sources 
with a luminosity exceeding 10 40 ergs/sec. Run 10 would 
therefore be a good candidate to explain the nature of the 
ULX in M82. 

In run 7, three stars experience tidal inspiral between 
3 < T < 4 Myrs and produce ULX sources before the stars 
are tidally disrupted by the IMBH or turned into compact 
remnants. At later times, the innermost stars are always too 
far away from the IMBH for significant mass transfer, hence 
run 7 is not able to explain the ULX in M82. In total we find 
that in 10 of the 12 performed runs a ULX source is produced 
at least once between 3 and 12 Myrs and that 4 out of the 
12 runs create an X-ray source brighter than 2- 10 39 erg/sec 
within the age range of MGG-11. The average time an X- 
ray source with a luminosity brighter than 10 39 erg/sec is 
present in our runs is 0.96 Myrs, while a source brighter than 
10 40 erg/sec is present for 0.62 Myrs. These values would 
be lower if tidal destruction of stars with super-Eddington 
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Figure 9. Same as Fig. [S] but now for run 7. This cluster also 
experiences several inspiral events between 3 and 12 Myrs, but 
none of them produces an X-ray source which is in the right age 
range and bright enough to explain M82 ULX-1. 

heating rates would be taken into account (see section Iq"31 . 
however, at the moment our code does not follow the RLOF 
phases which would significantly prolong the lifetime of the 
ULX sources and increase their luminosities. We therefore 
conclude that an IMBH in a dense star cluster like MGG-11 
has a high chance of creating an ULX and that our model 
can explain the ULX source in M82. 



6 CONCLUSIONS 

Observations of X-ray sources with luminosities higher than 
the Eddington luminosity of stellar mass objects provide 
strong support for the existence of IMBHs. However, the 
mere existence of an IMBH does not guarantee the exis- 
tence of ULXs; there has to be a mechanism which accounts 
for the accretion of gas by the IMBH. 

In this paper we confirm the hypothesis that massive 
(M* > 1000 Mq) stars, which may be the progenitors of 
IMBHs, form naturally in young dense clusters as a result 
of many collisions between young stars. We assume that 
the massive star indeed forms an IMBH, and we perform 
N-body simulations of the interaction between the IMBH 
and its host cluster, while accounting for tidal enc ounters 
between stars and the IMBH. Madhusudhan et al. f2005) 
found through stellar evolution studies that massive main- 
sequence stars with masses M* > 8Mq in orbits of 6 to 30 
times the radius of the donor star around an IMBH are nec- 
essary to explain ULX sources. Our simulations show that 
tidal heating of stars in young star clusters will create ex- 
actly such systems. We find that, as a result of tidal heating, 
stars are likely to circularise at distances of a few tidal radii 
from the IMBH. Once the star has circularised, accretion of 
stellar gas by the IMBH is sufficient to account for lumi- 
nosities as high as L x Si 10 39 erg/s. On average, the X-ray 
luminosities of the circularised stars exceed 10 39 ergs/sec for 
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almost 10 6 yrs within the first 12 Myrs in our runs, adding 
further credibility to the scenario that ULX are IMBHs ac- 
creting gas lost from nearby companion stars. 

Our conclusions differ from those reached by Blecha et 
The main reason could be that the IMBHs in our 
runs are more massive than the ones they considered and 
they found that higher-mass IMBHs have a higher chance 
of capturing companions. In addition, they also neglected 
the effect of mass segregation, whereas in our runs we find 
that mainly massive stars circularise around the IMBHs. 

As we noted in §©, stars which are tidally captured by 
an IMBH are typically among the most massive stars in the 
cluster. The massive star loses much of its energy to other 
stars by dynamical friction. In a single mass distribution 
(as is often assumed in an analytical treatment of the prob- 
lem^ diffusion in energy space is very slow (Bahcall & Wolf 
1976), but when a spectrum of masses is present, energy loss 
due to dynamical (non-elastic) interactions with field stars 
(as opposed to tidal interactions with the IMBH) plays an 
important role. This further emphasises the importance of 
mass segregation. 

In spite of the fact that some highly simplifying as- 
sumptions were made, the final estimate for the rate result- 
ing from our analytical treatment is of the same order of 
magnitude as what we find from our simulations. 

After the donor star leaves the main sequence, it forms 
a compact remnant. For high-mass donor stars the remnant 
will be a stellar mass BH and the subsequent evolution of the 
IMBH-BH binary is driven by the emission of gravitational 
waves until event ually th e binary members merge. Hopman 
& Portegies Zwart l(2005l) have shown that the event rate for 
this is likely to be high enough to be detectable by LISA, in 
particular if the IMBH mass is larger than ~ 3 x 10 3 Mq. 
Observations of gravitational waves from such a binary will 
give further support for the scenario discussed in this paper. 
In addition, if IMBHs are formed in star clusters near the 
centres of galaxies and spiral in due to dynamical friction 
on the field stars, mergers of the IMBHs with the MBHs 
at galactic centres w ould a lso be important LISA sources 
(Portegies Zwart et al l 2005F) . 

The mechanism of stellar capture by gravitational wave 
radiation is similar to that of tidal capture. Mass segregation 
probably plays a less significant role in the case of the inspi- 
ral due to gravitational wave emission by compact remnants 
when spiralling in to a MBH of M. > 1O 6 M0, because in 
that case the mass ratio of the field stars is less extreme than 
is the case for a young dense stellar cluster. The number of 
compact remnants in galactic centres is not known. Obser- 
vationally the mass of these objects can be constrained in 
our Galactic Centre by finding deviations from pure Kep- 
lerian orbits, in particular pericenter-shift (Mouawad et al. 
I 2004) . Current observational constraints are not conclusive, 
but this situation will improve with more accurate measure- 
ments of the stellar orbits in our Galactic Centre. 

In addition to direct disruption and tidal inspiral, we 
also find a large number of stars that experience one or sev- 
eral strong tidal encounters, but are scattered to wider orbits 
before they can circularise near the IMBH. As a result of the 
strong tidal encounters, such tidally heated stars or "tidally 
scattered" stars are expected to show signs of mixing, large 
spin and stripping; this may be directly observable in the 
Galactic centre (Alexander & Livi<E 2001). We found that on 



average there are 45 stars which have one or more interac- 
tions of less than 5rt with the MBH and are then scattered 
away. We expect that this number scales linearly with peri- 
center distance. 

We have assumed that the IMBHs in the clusters formed 
in a runaway merger of young stars, during a very early stage 
of cluster evolution. However, other scenarios of IMBH for- 
mation have been discussed in the literature. IMBHs may 
form a s remnants of population III stars (e.g. Madau & Rees 
I 20011) . or as the merger of stella r mass black holes due to 
GW emission (Miller & Hamilto rTioOl]) . If IMBHs indeed 
form in these scenarios, and they still reside in stellar clus- 
ters, the stars in these clusters will have a much lower max- 
imal mass than in the young clusters we discussed here. If 
a star is tidally captured, accretion will most likely not lead 
to high luminosity to account for the most luminous ULXs, 
although accretion via the subsequent RLOF will lead to a 
low-luminosity X-ray source. 
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